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Vortex trimer in three-component Bose-Einstein condensates 
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Vortex trimer is predicted in three-component Bose-Einstein condensates (BEC's) with internal 
coherent couplings. The molecule is made by three constituent vortices which are bounded by 
domain walls of the relative phases. We show that the shape and the size of the molecule can be 
controlled by changing the internal coherent couplings. 
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I. INTRODUCTION 



Recent advances in realizing Bose-Einstein conden- 
sates (BEC's) in ultracold atomic gases have opened new 
possibilities of quantum physics [l|, [2| . One of them is in- 
terpenetrating superfluids, a mixture of two or more su- 
perfhiids. Such multicomponcnt BEC's can be realized 
when more than one hyperfine spin state is simultane- 
ously populated or when more than one species of atoms 
are mixed. The s-wave scattering wave-length can be 
tuned via a Feshbach resonance 043 • Moreover, recent 
experimental achievement of a condensate ofytterbium 
offers condensations up to five components [6j. Stabil- 
ity condition of multicomponent BEC's was studied in 
Ref. [7|. One of the most important consequences of su- 
perfluidity is the existence of vortices. Vortices in mul- 
ticomponent BEC's have been realized experimentally 
[1, Hj], and structures of those vortices are much richer 
than those of single components [T34I3- 

In the case of multiple hyperfine spin states, the inter- 
nal coherent coupling between multiple components can 
be introduced by Rabi oscillations. This case is simi- 
lar to two gap superconductors with Josephson coupling 
between the two gaps. A sine- Gordon domain wall of a 
phase difference of two components is allowed [l3[ . More- 
over, an integer vortex is split into two fractional vortices 
with fractional circulations, and they are connected by a 
sine-Gordon domain wall with the total configuration be- 
ing a molecule of two constituent vortices, namely, a vor- 
tex dimer [l4| [l5| . Therefore it is natural to ask whether 
a molecule made of more than two vortices is possible in 
some case, or how domain walls connect among them if 
it is possible. 

In this paper we explicitly construct a vortex trimer, 
namely, a molecule made of three constituent vortices 
winding around respective three components of BEC's 
with internal coherent couplings induced by Rabi oscilla- 
tions. Varying the internal coherent couplings, the shape 
of the molecule is changed accordingly. We also find a 
dependence of the size of the vortex trimer on the mag- 
nitude of the Rabi frequency. 



II. THE GROSS-PITAEVSKII MODEL 

We consider three-component BEC's of atoms with 
equal mass m, described by the condensate wave func- 
tions tpi (i = 1,2,3) with the energy functional 



E 



E 



d 2 x 



(1) 



where atom-atom interactions are characterized by the 
coupling constants g.^ = gji, /x, is a chemical potential, 



and a symmetric tensor uji 



(uin — 0) stands for the 



Rabi frequency between the ith and jth components. In 
this paper, we consider the case with \i\ = ^ 2 = A*3 = Mi 
.9ii = .922 = .933 = 9, and g 12 = g 2 3 = 931 = 9 for 
simplicity, but the general case is straightforward. Wc 
also assume g + 2g > for the stability of ground states. 

In the following, we will separately study two cases: 
the case with g ^ g (det gij ^ 0) and the t/(3) symmetric 
case with g — g (det g^ = 0) . 

Let us first study the former case. When all the Rabi 
frequencies vanish, the ground state is given by 



Mi 



/i 



.9 + 2.9 



(i = 1,2,3). (2) 



The topology of the ground state is characterized 
7Ti[[/(l) 3 ] ='2fflZffil Once the small Rabi fre- 
quencies are turned on, only the overall U(l) symmetry 
remains contact and the homotopy group also reduces 
to 7Ti[Z7(l)] = Z. At the same time, the magnitudes 
of the condensates are modified. This is because the 
Rabi frequencies yield potentials on the relative phases 
of 8i = arg ipi . The ground state can be obtained by solv- 
ing a variational equation SE/Stpi — 0. We denote the 
condensate wave function of the ground state by 



il)i = Vie i6i , (Vi>0). 



(3) 



in Eq. ([T]) can be written as 
Oj). For uiij > 0, the phases Qi and 
9-i tend to coincide, Qi = Oj. to reduce this interaction 



The last term 
~2viVj0Jij cos(#i - 
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FIG. 1: (Color online) The parameter region where Eq. Q 
holds. The left panel shows the boundary surface on which 
Eq. (j4| holds inside (the region A) while it does not hold 
outside (the region B). The plot range of the left panel is 
—0.05 < oj-ij < 0.1. The right panel shows the cross sec- 
tion LJ23 = W31 . Some points { Ay , Bij , O} on the surface are 
shown: A12 is placed at (u)i2, 1^23, W31) = (—0.05, 0.1, 0.1), B12 
is at (0.1,0,0), and the other points are obtained by the Z3 
rotations around O. 
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energy. Therefore, when > hold for all i, j = 1, 2, 3, 
all the phases coincide, 

0i = o 2 = 9 3 , (4) 

in the ground state. We numerically find that relation 
(U]) is satisfied in the parameter region A shown in Fig. [T] 
We note that Eq. ((4]) holds not only when all wy's are 
positive but also in the region where one of the w^-'s 
is negative. In this paper, we consider region A, where 
relation (U) holds. The other region is frustrated, which 
will be studied elsewhere. 



FIG. 2: (Color online) The panels (a), (b), and (c) show the 
profiles of the density \ipi\ 2 , l^l 2 , and \1p3f f° r the (1,0,0)- 
vortex in the case where g 7^ g, respectively. The arrows 
show a phase vector (Re(|/>i),Im(^«)). The energy density is 
shown in the panel (d). We choose h = m = 1, g = 1000, g = 
900, /x = 100. The Rabi frequencies are CJ12 = UJ23 = 6J31 = 0. 
The sizes of the boundary are L = — 15 to L = 15 for both x 
and y. 



III. VORTEX CONFIGURATIONS 

The nontrivial first homotopy group immediately leads 
to the existence of superfluid vortices. In particular, 
the case with g 7^ g would have three different kinds 
of vortices because of the homotopy group 7Ti[i7(l) 3 ] = 
ZffiZffiZ (when Uij = 0). 

Let us consider an integer vortex configuration in 
which all the condensations ipi have unit winding in 
U(l)'s. The asymptotic behavior of such a configuration 
at large distance from the vortices should be 



(^1,^2,^3) -> (v 1 e ie ,v 2 e ie ,v 3 e^ 



(5) 



which satisfies the ground-state condition Here 9 

stands for the angular coordinate as x + iy = re 1 ® . 
We will show that this vortex is deformed to a vortex 
trimer made of three constituent vortices, (v\e lS , v 2 , U3), 
(vi, v 2 e ie , V3), and v 2l v 3 e t6 ), which we call (1,0,0)-, 
(0,1,0)- and (0,0,l)-vortices, respectively. Since it is a dy- 
namical problem if the constituent vortices make a bound 
state or not, let us see the two cases Wij = and Wjj 7^ 
separately. 



When the Rabi frequencies vanish = 0), the con- 
stituent vortices do not make a molecule since they re- 
pel each other. Instead, the constituent vortex can exist 



alone; see Fig. [21 where a numerical solution 18| of the 
constituent vortex is shown. The tension (energy per unit 

length) of the constituent vortex is given by — 
with system size L and healing length £. 



In I 



On the other hand, when the Rabi frequencies are 
not zero, the unit constituent vortex alone is unstable 
because semi- infinite domain walls are attached to it. 
This domain wall supplies attractive force between the 
constituent vortices and can be balanced with repulsion 
among them, so that the constituent vortices form a vor- 
tex trimer (vortex dimers are stable only when two of the 
Rabi frequencies are zero). To understand this better, it 
is useful to consider a reduced model from Eq. ((T|) by 
fixing the amplitudes ~ Vi (ui ~ v for simplicity). 
Then we are left with three phases: 



e 



2,3 



'2,3,1 
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(6) 
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The Hamiltonian of the reduced model is given by 



2„,2 



H 



6m 



(V9) 2 + (( v ^) 2 - w< cos Si) 



(7) 



where we have introduced the renormalized couplings 
£"1,2,3 = ^§2^23,31,12- This approximation is valid only 
when the Rabi frequencies are much smaller than the 
other coupling constants. [l9| For example, let us con- 
sider the (1, 0, 0)-vortex, with relative phases given by 



to=0, a a = -0i, to = 0i. 



Then the potential term reads 



V 



6m 



(<I>2 COS £ 2 + £3 COS to) ■ 



(8) 



(9) 



When 0)2,3 > 0, ^2 = 7r is an unstable point and a semi- 
infinite domain wall appears on the negative region of the 
real axis. Its tension is given by 



Ti 



??2 



T 2 

^31' 



T, 



8^ 



m(g + 2g) 



(10) 



This is the origin of the attractive force between the con- 
stituent vortices. Note that, since we have two relative 
phases to and to, one may naturally imagine two inde- 
pendent domain walls. Each domain wall has the tension 
T31 (when we set wi 2 = 0) and T12 (when we set W31 = 0). 
However, for the (l,0,0)-vortex, the two relative phases 
are related as to — —to an d these two domain walls stick 
together and form a bound state. Indeed, the total ten- 
sion 7i is the square root of the sum of T 2 2 and T| 2 as 
shown in Eq. (TIT)]) , which is smaller than the sum of the 
two tensions, T x < T 12 + T 31 . 



IV. NUMERICAL SOLUTIONS OF VORTEX 
TRIMERS 

We have numerically found stable vortex trimers as 
unique solutions under the boundary condition ([5]) in a 
wide range of the parameter region A in Fig. [TJ Here in 
Fig. [3] we show several numerical solutions as examples 
with g = 1000 and g = 900 (/z = 100 and m = h = 1). 
The rightmost panels show the energy density in which 
the partonic structure is clearly seen. The contours 
therein show contributions from the last term of Eq. ((T|) . 
Since the distance between the constituent vortices are 
close, we cannot see domain walls. Nevertheless, quali- 
tative estimation from the reduced model is quite useful, 
as will be seen below. 

Each line of Fig. [3] gives a molecule with different Rabi 
frequencies. First of all, the fourth line of Fig. [3] shows 
a Z3 symmetric trimer where the Rabi frequencies are 
all equal as wi 2 = W23 = ^31 = 0.05 . One can see 
that the phases at the spatial infinity are indeed aligned 
(0i = 02 = 03)- Next, by changing wi 2 from the sym- 
metric case, we can observe how the shape of the trimer 



is deformed. Since the Rabi frequency uji 2 controls the 
interaction between the (1,0,0)- and the (0,l,0)-vortices, 
the equilateral triangle is deformed to an isosceles trian- 
gle. The third line of Fig. [3] shows the vortex trimers 
with UJ12 = 0.01, in which the attractive force between 
the (1,0,0)- and the (0,l,0)-vortices is smaller than those 
between the other two pairs. Therefore the internal an- 
gle at the vertex at the (0,0,l)-vortex is larger than 7r/3. 
We also show the vortex trimer when W12 = in the 
second line of Fig. [3] Since no attractive force exists be- 
tween the (1,0,0)- and the (0,l,0)-vortices, the shape of 
the molecule becomes a stick as expected. We also find 
the molecule even when cji 2 is negative (= —0.01) while 
W23 = W31 = 0.05; see the first line of Fig. [3] We still ob- 
serve a stick-type molecule whose length is slightly larger 
than that for wi 2 = 0. In the last two lines of Fig. [3J we 
have chosen UJ12 = 0.2 and 0.5, which are larger than 
CJ23 = W31 = 0.05 (the Z3-symmetric case). Since the at- 
tractive force between the (1,0,0)- and the (0,l,0)-vortices 
is stronger than those for the other two pairs, we see that 
the corresponding edge of the triangle becomes shorter 
than the other two edges. Since wi 2 = 0.5 yields an at- 
tractive force that is too strong, the triangle collapses, as 
shown in the last line of Fig. [3j 

We have seen that the shape of the triangle changes 
according to the choice of the Rabi frequencies. Here, 
we investigate a correlation of the Rabi frequencies, set 
to be equal uji 2 = W23 = ^31 = and the size of the 
equilateral triangle, see Fig. 2) We numerically find the 
relation 



A ~ 0.56 oj- 



(11) 



for the range 0.01 < ui < 0.1, where A stands for the 
length of the edge of the equilateral triangle. 

Let us finally consider a vortex trimer in the C/(3) sym- 
metric case (g = g) where all of the terms in Eq. (TTJ) ex- 
cept for the last one are invariant under ip — ¥ U?p with 
U E U(3). A striking difference from the previous case 



with g ^ g can be best seen in the limit where uii 



0. 



The ground state is degenerate and its order parameter 
space is U{3)/U(2) ~ S 5 defined by £? =1 \ip t \ 2 = The 
first homotopy group of the ground state is trivial, and 
there are no topologically stable vortices. When the Rabi 
frequencies are not zero, the order parameter space be- 
comes U(l) implying the existence of stable vortex con- 
figuration. Fig. [5] shows numerical solutions for several 
choice of the Rabi frequencies. We examine two choices 
for the Rabi frequencies with g = g = 1000 (/x = 100 
and h = m = 1): i) u)\2 = lu 2 3 = LJ31 = 0.1, which 
leads to the ground state condensation (\ipi\, IV^I, 1^3 1) = 
(0.183,0.183,0.183), and ii) w 2 3 = 0.05, lj 12 = w 3 i = 0.1, 
which leads to (|^i|, \i/j 2 \, \ipa\) = (0.203,0.171,0.171). 
Although the profiles of the condensation are not ax- 
isymmetric, the total energy density is universally ax- 
isymmctric; see the rightmost panels in Fig. [5] Unlike 
the previous case with g 7^ g, it is impossible to see a 
partonic nature from the energy density when g = g. 
These things are related to the fact that there are no con- 
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FIG. 3: (Color online) The left three panels show the profiles of the density \ipi\ 2 and the phases (Re[0i], Im[^]i) for the unit 
vortex trimer in the case g ^ g, respectively. The rightmost panel shows energy density and the contour corresponds to the 
Rabi potential. The constants are taken as h — m = 1, [i = 100, g — 1000, g — 900, and ct/23 = ^31 = 0.05. We change the 
Rabi frequency a/12 from the top to the bottom as a/12 = —0.01,0,0.01,0.05,0.2,0.5, respectively. The sizes of the boundary 
are L = — 8 to L = 8 for both x and y. 
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FIG. 4: (Color online) Loglog-plot of |wi2| vs the length (^4) 
of an edge of the equilateral triangle with uj\2 = ^23 = U31 ■ 
The corresponding configuration is given in the fourth line of 
Fig. [3] The parameters are fixed as h = m — 1, fj, = 100 and 
g = 1000 and g = 900. 
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FIG. 5: (Color online) The panels (a), (b), and (c) show the 
profiles of the density l^il 2 ) | "02 1 2 , and \ip3\ 2 for the minimal 
vortex in the !7(3)-symmetric case, respectively. The energy 
density is shown in panel (d). We choose ft = m = l,g — 
g = 1000, /i = 100. The Rabi frequencies are chosen as CJ12 = 
[J23 = W31 = 0.1 in the first row, and UJ23 ~ 0.05, 0^12 = UJ31 
1 in the second row. The sizes of the boundary are L — — 15 
to L — 15 for both x and y. 



stituent vortices standing alone in the limit — > 0. This 
configuration can be regarded as a skyrmion in the CP 2 
nonlinear sigma model with the two-dimensional complex 



projective space CP 2 ~ S^S 1 ~ SU(3)/[SU(2) x [7(1)], 
instead of CP 1 ~ S 3 /S 1 ~ SU(2)/U(1) ~ S 2 for two 
component BECs [H |. 



V. DISCUSSION 

Finally, we comment on the possibility of realization in 
experiments. In this paper, we have solved the GP equa- 
tion in the static system by fixing the phase winding and 
the constant density at the boundaries [18]. This implies 
that the vortex trimer is created in rotating BEC's in lab- 
oratory experiments. First, an integer vortex is created 
as usual by gradually increasing the rotation, and then 
it will be split into the vortex trimer. Two component 
BEC's of different hyperfine states of the same atom have 
been already realized using the |1, — l) and |2, l) states 
H| and the |2, l) and |2,2) states [l6( of 87 Rb, respec- 
tively. Three component system should be possible using 
a mixture of those states of 87 Rb. For that, one needs to 
use the optical trap which has been recently been realized 
in a two-component system [l7| . 

In this paper, we have investigated the vortex trimer 
in the three-component BEC's. We expect that the vor- 
tex A^-omers can also be constructed in ./V-component 
BEC's. However, for N > 4, the choice of the internal 
coherent couplings becomes nontrivial. For instance, 
for N — 4, we expect that the choice UJ12 = UJ23 = W34 = 
UJ41 > gives a symmetric tetramer, but we still have 
additional parameters W13 and W24. We will study the 
7V-omer elsewhere. 
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